KU-TP 057 



Charged Black Holes in String Theory with Gauss-Bonnet Correction in Various 

Dimensions 

Nobuyoshi Ohta 1 '* and Takashi Torii 2 ' ^ 

department of Physics, Kinki University, Higashi-Osaka, Osaka 577-8502, Japan 
^Department of General Education, Osaka Institute of Technology, Asahi-ku, Osaka 535-8585, Japan 

(Dated: October 24, 2012) 

We study charged black hole solutions in Einstein-Gauss-Bonnet theory with the dilaton field 
which is the low-energy effective theory of the heterotic string. The spacetime is D-dimensional and 
assumed to be static and spherically symmetric with the (D — 2)-dimensional constant curvature 
space and asymptotically flat. The system of the basic equations is complex and the solutions are 
obtained numerically. We identify the allowed parameter region where the black hole solutions exist, 
and show configurations of the field functions in D — 4 - 6 and 10. We also show the relations of the 
physical quantities of the black holes such as the horizon radius, the mass, the temperature, and so 
^ ' on, and find several results. The forms of the allowed parameter regions are different depending on 

, the dimension. There is no extreme black hole solution with T = that can be obtained by taking 

the limit of the non-extreme solutions within the parameter range we chose. Entropy of the black 
' holes in the dilatonic theory is always larger than that in the non-dilatonic theory. Our analysis 

includes the higher order term of the dilaton field which is not in our previous works. Its effect 
remarkably appears in five dimensions and is given in the appendix. 
CO 1 By our analysis it is found that the properties of the black hole solutions strongly depend on the 

dimension, charge, existence of the dilaton field. Hence both the detailed analyses of the individual 
systems and the investigations from the systematic point of view are important. 
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I. INTRODUCTION 



One of the most important problems in theoretical physics is to find the quantum theory of gravity and to apply it 
to physical systems to understand physics at strong gravity such as a black hole singularity. The leading candidates 
. for that including all the fundamental forces of elementary particles are the ten-dimensional superstring theories or 
\Q ■ eleven-dimensional M-theory. The area where such quantum gravity plays the significant role includes the cosmology 
and black hole physics. There has been interest in applications of string/M theories to these subjects. Since it is still 
difficult to study geometrical settings in string/M theories directly, most analyses have been made by using low-energy 
^vq | effective theories inspired by string/M theories. The effective theories are the supcrgravities which typically involve 
not only the metric but also the dilaton field (as well as several gauge fields). 

The first attempt at understanding black holes in the Einstein-Maxwell-dilaton system was made in Refs. [1, 2], in 
which a static spherically symmetric black hole solution with a dilaton hair was found. After that, many solutions 
were discussed in various models. On the other hand, it is known that there are correction terms of higher orders in 
the curvature to the lowest effective supergravity action coming from the superstring theories. The simplest correction 
is the Gauss-Bonnet (GB) term coupled to the dilaton field in the low-energy effective heterotic string [3]. It is then 
natural to ask how the black hole solutions arc affected by the higher order terms in the effective theories. 

When the dilaton field is dropped or is set to a constant by hand, the total action consists of the cosmological 
constant, Einstein- Hilbert and GB term, which are the first three terms in the Lovelock theory [4, 5]. Motivated 
by this observation, there have been many works on the black hole solutions in the Lovelock theory [6-9]. In the 
four-dimensional spacetime, the GB term does not give any contribution because it becomes a surface term and gives 
a topological invariant. Boulware and Deser (BD) [10] discovered a static, spherically symmetric black hole solutions 
of such models in more than four dimensions. In the system with a negative cosmological constant, black holes can 
have horizons with nonspherical topology such as torus, hyperboloid, and other compactificd submanifolds. These 
solutions were originally found in general relativity and are called topological black holes [11]. Topological black hole 
solutions were studied also in the Einstein-GB (EGB) theory [12]. It is very interesting to see how these arc modified 
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by the presence of a dilaton field. One of the purposes of this paper is to study charged black hole solutions with 
higher order corrections as well as the dilaton field. 

There has been recently a renewed interest in these solutions as an application to the calculation of shear viscosity 
in the strongly coupled gauge theories using the black hole solutions in the five-dimensional EGB theory via AdS/CFT 
correspondence [13, 14]. Almost all these studies considered the pure GB term without the dilaton field, or assumed 
a constant dilaton, which is not a solution of the heterotic string theory. It is, however, expected that AdS/CFT 
correspondence is valid within the effective theories of supcrstring which necessarily involve the dilaton field. The 
inclusion of the dilaton field was also considered by Boulware and Dcser [10], but exact black hole solutions and their 
thermodynamic properties were not discussed. Callan et al. [15] considered black hole solutions in the theory with a 
higher-curvature term R pvpa R^ vpcr and the dilaton field, and Refs. [8, 9] took both the GB term and the dilaton field 
into account in four-dimensional spacetime. 

Hence it is important to study black hole solutions and their properties in the theory with both the higher order 
corrections and the dilaton field in general dimensions. In our earlier paper [16], we focused on asymptotically fiat 
solutions and studied a system with the GB correction term and the dilaton field without the cosmological constant 
in various dimensions from 4 to 10. They are spherically symmetric with the (D — 2)-dimensional hypersurface of 
curvature signature k — 1. We have then turned to planer symmetric (k — 0) topological black holes, but have found 
that no solution exists without the cosmological constant. In the string perspective, it is more interesting to examine 
asymptotically anti-de Sitter (AdS) black hole solutions with possible application to the AdS/CFT correspondence in 
mind. So in the sequel paper [17], we have presented the results on black hole solutions with a negative cosmological 
constant with k = 0. In fact, shear viscosity has been computed using the result in that paper [14]. Other cases 
of solutions with k = ±1 and a negative cosmological term as well as topological black holes with and without the 
cosmological constant were studied in [18, 19], and their global structures were studied in Ref. [20]. Cosmological 
solutions were also considered in Ref. [21]. Other solutions are discussed in Refs. [22-29]. 

In the study of superconductors and superfluidity using the AdS/CFT correspondence, charged black hole solutions 
play important roles [30]. Hence it is also significant to extend our above studies of the neutral black hole solutions 
to charged ones. In particular, we should consider the inclusion of the dilaton field again. 

There is another reason for our study of the system. In our above study of the neutral black holes [16-20], we did 
not consider the higher derivative term of the dilaton field, which also appears naturally as we will see below [3] . We 
intend to incorporate this term also and study how this modifies the solutions. As a result of our analysis, however, 
there is not much qualitative difference. Thus the results in our earlier papers should be useful. 

This paper is organized as follows. In Sec. II, we first present the action and give basic equations to solve for the 
system of Einstein-Maxwcll-Gauss-Bonnct term coupled to the dilaton field with a cosmological constant, although we 
will focus on the asymptotically flat solutions without cosmological constant later on. The case with the cosmological 
constant is left for future study. Boundary conditions and symmetry properties of the theory are also discussed in 
order to apply them in our following analysis. In Sec. Ill, we briefly summarize the thermodynamic properties of the 
solutions in the theory. In Sec. IV, we also briefly review the results for non-dilatonic solutions to see the difference 
from the dilatonic ones. We then proceed to the study of the charged black holes for D = 4, 5, 6 and 10 in Sees. V, 
VI, VII and VIII, respectively. The study of the effects of the higher derivative of the dilaton field for the neutral 
black holes is delegated to the appendix, because we find that there is not much difference in the properties of the 
black hole solutions except for five dimensions. Some detailed discussions are given there. Section IX is devoted to 
conclusions and discussions. 



II. DILATONIC EINSTEIN-MAXWELL-GAUSS-BONNET THEORY 



A. Action and basic equations 



We consider the following low-energy effective action for the heterotic string theory in one scheme [3]: 
1 
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where «f is a D-dimcnsional gravitational constant, cf> a dilaton field, F a gauge field strength, H a 3-form, a 2 = o 1 



no 

is a numerical coefficient given in terms of the Regge slope parameter a', and i? GB 
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Hfj, pa H vpa . In the original derivation of the effective action, it was first 



is the GB term. H 2 

derived in the Einstein frame from the S-matrix calculation in the string theory, and then transformed into the string 
frame [3]. It is common and convenient to interpret results in the Einstein frame. Hence we transform Eq. (2.1) into 
the Einstein frame, reduce to D dimensions, and use the field redefinition ambiguity [3, 31] 



5<f> = a , [b 1 R + b 2 {d4>) 2 + b 3 V 2 (t>}, 
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where 7 = 1/2. We have set H = because we focus on the effects of the dilaton field and the gauge field. Note that 
H = is a solution of the field equations. 

When the action (2.4) is derived from (2.1), a numerical constant a in the action (2.4) should be a = 1. In our 
previous papers [16-20], however, we have examined the system with a = 0, which corresponds to the case without 
the higher order term of the dilaton field. To see how the solutions are different with and without the higher order 
term of the dilaton field, we have included the constant a, and here we study mainly a = 1 case. The difference 
between a = and a = 1 is discussed in the Appendix. We have also included a cosmological constant with the 
dilaton coupling A although we study asymptotically flat solutions in this paper. 

In the above process, higher order terms (> 0(a' 2 )) are dropped. This is allowed because the effective low-energy 
action can be determined up to the field redefinition when it is read off from the scattering amplitudes computed 
in the string theories. Also note that the effective action was originally computed in the Einstein frame. All this 
means that there is no absolutely preferred form of the action if they are related up to terms of order a' 2 . There 
might be some quantitative differences if we adopt theory of different choice of the coefficients, but we expect that 
the qualitative properties will not change significantly. It is the system (2.4) with a = 1 that we study in this paper. 

Let us consider the metric and field strength 
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where hijdx l dx J represents the line element of a (D — 2)-dimcnsional constant curvature space with curvature (D 
2){D - 3)fc and volume E fe for k = ±1, 0. 

The field equations following from Eq. (2.4) are 
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where we have defined the dimcnsionlcss variables: f = rf-y/a^, A = o^A, and the primes in the field equations denote 
the derivatives with respect to r. Namely we measure our length in the unit of v / 02- In what follows, we omit tilde 
on the variables for simplicity. We have also defined 
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The field equation for the Maxwell field (2.9) is easily integrated to give 
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where q is a constant corresponding to the charge. It should be noted that the quantity q in this equation is actually 
the rescaled one q = q/a( D ~ 3 ^ 2 related to the real charge q. But as we mentioned before, tilde is omitted here and 
below. 

Our task is reduced to setting boundary conditions for the metric functions B, 5 and the dilaton field <p and integrate 
the above set of equations, just like in our previous papers [16-20]. 



B. Boundary conditions and asymptotic behaviors 



In this paper we consider non-extreme solutions which has a non-degenerate black hole horizon. Hence at the 
horizon r# , 



B H = 0, B' h ^0, 
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where Bh = B(th)- Here and in what follows, quantities evaluated at the horizon are represented by a subscript H. 
At the horizon, Eq. (2.6) gives 
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Henceforth we restrict our consideration to k = 1 and A = case. Combining Eq. (2.16) with Eq. (2.8) evaluated at 
the horizon, we obtain the quadratic equation determining cj)' H : 

A 2 $ 2 + A 1 <j>' H + A = (2.18) 

where 

A 2 = C^r 2 I P +2 [-2{D-3)lc((D-2)^ 2 [2C{D-3)-l]+D-A^ +l}q 2 r 4 H 

+C(D - 2)rjf{c{D - 4)(C(D - 2) 7 2 [2C(L> - 5)(D - 3) + 3D - 11] + C(D - 4) 5 + 3D - ll) + 2(D - 3)} 

(2.19) 
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+2(L> - 3)g 2 r|f +5 {2C , (C - 1)(D - 2) 7 2 [2C(L> - 4) + 1] + [C(D - 4) + l] 2 } + 4C(D - 3) 2 q 4 r 9 H j 2 (2.20) 

A = -\C 2 {D - 2) 2 (D - lh{c(D - A)[C{D - 4)(D + 1) + 4] - 2(D - 2)}r|f + 2(D - 3) 2 q 4 r 8 H j 

-2(D - 2) 3 q 2 1 [3(D - 4)C 2 + 6C - l]r 2 ° +i (2.21) 

The (d4>) 4 term in Eq. (2.4) does not contribute to the boundary condition of the dilaton field at the horizon. 
In the asymptotic region of r — > oo, we assume 

2M 

B^1-^L + ... (2.22) 

(5^0 (2.23) 
0^0 (2.24) 

where M is a constant corresponding to the mass of the black hole. Although 6 — > 600 and <f> — > 0oo, where 60a and 
4>oo are constant, these constants can be rescaled out to satisfy Eqs. (2.23) and (2.24) by the symmetries shown below, 
and the solutions are asymptotically flat. 

C. Symmetry and scaling 



It is useful to consider several symmetries of our field equations (or our model). 
First, our field equations (2.6) - (2.8) have a shift symmetry: 



<j> tf> - 0*, r^e^-/ 2 r, q -> e (D~^<f,/2 q ^ (A — s- e^ -7 -^* A), (2.25) 

where 0* is an arbitrary constant. 1 This changes the magnitude of the cosmological constant when we consider black 
hole solutions in its presence. Hence this may be used to generate solutions for different values of the cosmological 
constant, given a solution with some value of cosmological constant. Even without the cosmological constant, this 
symmetry can be used to change the asymptotic value of the dilaton field. 
The second one is another shift symmetry under 

d^>5-6*, t^e~ s 't, (2.26) 

with an arbitrary constant <5*, which may be used to shift the asymptotic value of S to zero. 2 



1 There are typos in our previous papers [18, 19] in the exponent of the transformation rule in the second term; they should have the + 
sign. 

2 For k = 0, there is another symmetry. The field equations (2.6)-(2.8) are invariant under the scaling transformation B — > c 2 B, r — > cr, 
with an arbitrary constant c. If a black hole solution with the horizon radius ru is obtained, we can generate solutions with different 
horizon radii but the same A by this scaling transformation. 
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III. THERMODYNAMIC AL VARIABLES 

Here we briefly summarize thermodynamical quantities of black holes to be used in the following discussions. The 
Hawking temperature is given by the periodicity of the Euclidean time on the horizon as (keeping k arbitrary) 
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Along the definition of entropy in Ref. [32], which originates from the Noether charge associated with the diffco- 
morphism invariance of the system, we obtain 



-2tt 



dC 



dR 



(3.2) 



where £ is the event horizon (D — 2)-surface, C is the Lagrangian density, e^ v denotes the volume element binormal 
to S. This entropy has desirable properties such that it obeys the first law of black hole thermodynamics and that it 
is expected to obey even the second law [33]. For our present model, this gives 
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It is noted again that is the volume of the unit constant curvature space, and £1 = y({d-i)/2) ' 
make the entropy non-negative [34]. 



(3.3) 



Smin is added to 



IV. NON-DIL ATONIC BLACK HOLE SOLUTION 



It will be instructive to compare our results with the non-dilatonic case. When the dilaton field is absent (i.e., 
Einstein-Maxwell-GB system), we substitute = and 7 = into Eqs. (2.6) and (2.7). In the D = 4 case, the GB 
term is total divergence and does not give any contribution to the field equations. As a result, the solution reduces 
to the Reissner-Nordstrom (RN) solution. 

For D > 5, the field equations can be integrated to yield [22] 
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and M is an integration constant corresponding to the asymptotic value m(oo) for the plus sign in Eq. (4.3). In the 
«2 — > limit, the solutions with the plus sign approach the RN solutions. This means that they can be considered to 
be the solution with GB correction to general relativity (GR). On the other hand, the solutions with the minus sign 
do not have such a limit. For these reasons, we call the solutions with plus (minus) sign the (non-)GR branch. 
For M = and 9 = 0, the metric function becomes 



B 




(GR branch) 
(non-GR branch) 



(4.4) 



where £ 2 S = (D — 3)4. Hence the spacetime is Minkowski in the GR branch while the spacetime is anti-de Sitter in 
the non-GR branch although the cosmological constant A is absent. 
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For the charged solution, besides the central singularity at r = 0, there can be another known as the branch 
singularity at finite radius rj, > 0, which is obtained by the condition that the inside of the square root in Eq. (4.3) 
vanishes. We find the M-r^ relation 



-D-l 



M 



8(£>-3) 4 



(D - 4)g 2 



8(D - 2)r 2(D ~ 2) 



(4.5) 



This implies that the branch singularity can appear for positive mass in charged black hole, but only for negative 
mass parameter in the neutral black hole. 

It can be shown that there is no black hole solution in the non-GR branch. On the other hand, in the GR branch, 
Eq. (4.3) evaluated at the horizon B = gives 
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8(D - 2) 3 r% D 
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This is the M-rn relation for the black hole without the dilaton field. For q = 0, we see that D = 5 is special because 
only in this case the mass goes to a finite value in the limit of r# — > 0, whereas it vanishes for other dimensions 
including four. This limit does not exist for charged ones (q =/= 0). 



V. D = 4 BLACK HOLE 

For D = 4, Eq. (2.18) to determine 4>' H reduces to 

Cr% 7 {2Cr% - q 2 [2C(2C - 1) 7 2 + 1] }<t>H + r H {2Cq 4 1 2 + q 2 r% [4C(C - 1) 7 2 + l] - lCr%)<t>' H 

+7 [l2C 2 r% - 2(6C - l)q 2 r 4 H + q 4 ] = 0. (5.1) 

The discriminant for 7 = i is 

Kl (q 2 - 2Cr%) 2 [C 2 q 4 + 2C(6C 2 - 3C + A)q 2 r A H + 4(1 - 6C 2 )r s H ] . (5.2) 

Our procedure for obtaining solutions is as follows. First we choose the suitable values of parameter q and the 
boundary condition (j>H (say 4>h = 0) at the horizon. Given the horizon radius r# together with (j)' H determined by 
(5.1), we integrate the field equations (2.6) - (2.8) outward from the horizon numerically There are two possible 
solutions for cj>' H in Eq. (5.1), but only the smaller one gives the black hole solutions. Although the dilaton field takes 
nonvanishing value in the asymptotic region gencrically, it can be set to zero by using the shift symmetry (2.25). 
Note that the variables defined by 

x = e*+«' 2 rH, y = eP- 2 ^' 2 q, (5.3) 

do not change under this shift symmetry. Then we redo this procedure by changing rjy but keeping the values of 
charge q and tfiu unchanged, which means that y is also unchanged, and obtain another solution with the different 
asymptotic value of the dilaton field (f)^. Using the shift symmetry to set (j)^ to zero, we have the solution with the 
same y but different x. In this way we obtain various solutions for one chosen y. Next, we repeat this procedure for 
different y (practically, with different q and the same 4>h)- After this process, we not only obtain solutions but also 
identify the parameter ranges where the solutions exist on the x-y plane. 

Since the values of x and y are unchanged by the shift symmetry (2.25), it is convenient to draw various diagrams 
in terms of these variables, x and y will be referred to as the "scaled horizon radius" and the "scaled charge", 
respectively. After identifying the parameter region where solutions exist, we can change these "scaled variables" to 
physical ones. 

We find that there are black hole solutions in the region shadowed by horizontal thin (blue) lines in Fig. 1(a). The 
horizontal and vertical axes are x and y, respectively. These horizontal (blue) lines actually consist of the sequences 
of solutions obtained by the above procedure. Although the lines are drawn for separate distance in y with intervals 
Ay = 0.2, it should be understood that there are solutions in the whole shadowed region. We call this region (where 
the black hole solutions exist) the allowed region. On the dashed (red) curve F determined by 



F: y 2 = -2x A + 3x 2 - 12 + ^3(x 8 - Ax e + 27x 4 - 24x 2 + 48) x 2 , (5.4) 
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1.5 2 2.5 3 2 2.5 3 3.5 4 



(c) (d) 
FIG. 1: The parameter region where the black hole solution exists in D — 4. The black hole solution exists in the allowed region 
which is shadowed by (blue) thin lines, (a) The horizontal axis is x = e 7 * H ' 2 rH and the vertical axis is y = e~ , ^ H q. On the 
dashed (red) curve F, the discriminant (5.2) is zero. Below this line, the discriminant is negative and there is no appropriate 
boundary condition on the event horizon. On the curve S the first derivative of the dilaton field on the horizon <j>' H diverges 
and the solution becomes singular. The line E represents the extreme solution. On the curve Bi, (f)^ diverges numerically 
which means (f>u diverges physically by the shift symmetry, (b) The magnified diagram around the top of the forbidden region. 
Around the curve B2 numerical calculation becomes unstable and we cannot find the black hole solution above it. On the curve 
B3, the third derivative of the dilaton field diverges while the Kretschmann invariant does not. (c) The allowed region in x and 
the physical temperature T of the black hole. The temperature becomes zero on the line Bi . The reason is that the horizon 
of the solutions on it have infinite size. The solutions on the boundaries F, B2, and B3 have non-zero finite temperature, (d) 
The diagram of the allowed region in terms of the physical quantity th and q. For each value of the charge q, there is a lower 
bound for the horizon radius th ■ 



the discriminant (5.2) vanishes. Below the curve F, the discriminant is negative and the values of (£>' H becomes 
imaginary, and hence there is no black hole solution. Furthermore, on the curve F, it can be found analytically that 
the second derivative of the dilaton field diverges at the horizon. The Kretschmann invariant 

= [B" - 3B'S' + 2B(5' 2 - 5")} 2 + 2(1 ^ ] (B' 2 - 2BB'S' + 2B 2 5' 2 ) + 2(£> ~ 2)3 (fc - B) 2 , (5.5) 

also diverges there, and the solution is singular. We call the region on and below the curve F the forbidden region. 
For q = 0, the boundary is x = 24 1 / 4 tn 2.213. For y > 3.0523 (above the top of the curve F), the discriminant is 
positive for all x (and rjr). 

The extremal condition given by Tjj = (or B' H = 0) is written in general dimensions as 

E: y 2 = 2[(D-2) 3 x 2 (°-V + (D - 2) 5 x 2 ^}. (5.6) 

In four dimensions, this reduces to 



E : y = 2x, 



(5.7) 
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and is depicted by the dotted (green) line E. Above the line E, B' H is negative and the horizon m is not a black hole 
horizon (i.e., a closed trapped surface). Actually we find that there is no black hole solution beyond the solid (black) 
curve S, on which <fi' H diverges. Hence the curve S is expected to give a boundary of the allowed region (for y < 4). 
The equation of S can be read off from Eq. (5.1), as the condition that the coefficient of <j>'^ vanishes: 



For small x { < 1), this singular curve S almost overlaps with the curve F. And it crosses with the extreme curve E at 



For y > 4, we find that the black hole solutions exist to the right of the boundary Bi. On Bi, we find that the 
asymptotic value of the dilaton field 4>oo diverges in the numerical integration, which means that 4>h diverges after <j>u 
is shifted such that 4>oa — > by the shift symmetry. Figure 1(b) is a magnified diagram around the top of forbidden 
region. For 3.0523 < y < 4, we find that there seems to be a boundary B2 slightly to the right of the curve S. We 
have confirmed that the black hole solutions exist on the right side of B2. In our numerical analysis, the calculation 
becomes unstable and stops for parameters just out of the event horizon before we let the parameters reach the curve 
S, which strongly suggests that there exists the boundary B2. On the other hand, expanding the field equations 
around r#, we do not find any singular behavior in the right region of S including B2 analytically. Unfortunately 
within our numerical accuracy, it is difficult to determine precisely where the true boundary is. However, there is 
certainly the boundary B2 different from S around y ~ 2.8, and it is plausible that it merges with S at the point 



For y < 3.0523, the allowed region splits into two parts. In the region to the right of the curve F, the black hole 
solutions exist. On the left side of the curve F, we also have what we call solitary solutions. We find that the right 
boundary of the region where the solitary solutions exist is partly F and partly B3 (see Fig. 1(b)), where the third 
derivative of the dilaton field diverges at certain radius r > m while the Kretschmann invariant JC does not. The 
endpoint of B2 and B3 is depicted by a (black) dot. Below this, there are black hole solutions only to the right of F. 

In Fig. 1(c), we show the temperature vs. the scaled horizon radius relation. There are solutions in the region 
where (blue) thin curves are plotted. It should be noted that the temperature vanishes on the boundary Bi, although 
Bi is different from the extreme curve E. The reason why the temperature vanishes is that the asymptotic value of 
the dilaton field 4>oo gets infinite and the horizon radius rjj and the charge q of the black hole on Bi become infinitely 
large after using the shift symmetry. The physical value of 4>h diverges and the solutions on Bi is singular. On the 
left lower boundary B2 around x = 1.6 ~ 2 the temperature is non-zero. The right-upper boundary is determined by 
q = 0, beyond which there is no solution. 

The diagrams drawn in terms of x and y are useful to identify the parameter region numerically where the black 
hole solutions exist. To discuss the physical properties, however, it is more convenient to show the diagrams in terms 
of the physical quantities such as rjj and q directly. In Fig. 1(d), we display the allowed region in terms of the 
charge q and the horizon radius r#. There is a lower bound on the horizon radius for any value of charge q. When 
the charge is small (g<5), the curve F gives the boundary. We can see that the minimum size of the black hole is 
almost independent of the charge in this range. When the charge is in the range 5<g<6.3, the curve F is again the 
boundary, but it is not the lower bound; there are solitary solutions in the thin region in the left part of the diagram. 
This cusp structure is surrounded by B2 and B3. The structure is so thin that it is difficult to create the black hole 
in this region in physical processes. For q > 6.3, the lower boundary is given by B2. The curve Bi with y > 2 in 
Fig. 1(a) is pushed to the upper right infinitely (th — > 00, q — > 00). 

So far we have presented the discussion of the parameter region where black hole solutions exist. We now focus 
on how the field functions of the black hole solutions behave for a typical boundary condition. Figure 2 shows the 
behaviors of the function <j>, m and 5 for the horizon radius ru = 3.2 for neutral and charged cases with q = 0, 5 
and 10. This horizon radius is chosen such that the solutions exist for these charges. For q = 10, the parameters are 
close to the boundary B2 in the diagram Fig. 1(d). We see that these functions have smooth behaviors. The dilaton 
field 4> increases monotonically for each charge (Fig. 2(a)). The mass function m of the neutral solution decreases 
near the horizon and increases towards a finite value as r increases [16]. On the other hand the mass function of 
the charged solution increases all the range from the horizon (Fig. 2(b)). The lapse function 5 decreases for small q. 
As the charge becomes large and the parameters approach B2, S increases around the horizon and decreases to zero. 
Since the dilaton field couples to the gauge field and the GB term, the effects of the dilaton field and higher curvature 
terms become significant as the charge becomes large. 

Up to this point, our discussions are rather qualitative with only figures of regions and behaviors of functions. To 
evaluate actually some quantities, it is necessary to have quantitative results. In order to get some idea on what are 
the typical physical quantities, here we tabulate them for the black hole solutions with the charge q = 0, 5 and 10 in 




(5.8) 



(x,y) = (2,4). 



(x, y) = (2, 4). 



Table I. 
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FIG. 2: Configurations of (a) the dilaton field 0, (b) the mass function m, and (c) the lapse function S of the black hole solutions 
in D — 4. The horizon radii are ru = 3.2, and the charges are q — (solid (black) line), q = 5 (dashed (red) line), q — 10 
(dotted (blue) line). 



q 




M 


Sh 


4>H 


<f>H 


T 







2.8 


1.50205 


0.291474 


-0.650704 


0.504050 


0.0282842 


3.34452 


3.2 


1.65494 


0.104926 


-0.439547 


0.253039 


0.0248377 


3.80579 


4.0 


2.02371 


0.0319963 


-0.256190 


0.110000 


0.0198896 


5.13666 


6.0 


3.00621 


0.00521331 


-0.106409 


0.0294481 


0.0132626 


10.0546 


8.0 


4.00251 


0.00156277 


-0.0586656 


0.0120864 


0.00994713 


17.0298 


5 


2.8 


2.48068 


0.502017 


-1.10819 


0.564272 


0.0143593 


3.70037 


3.2 


2.55368 


0.234784 


-0.871222 


0.347580 


0.0143019 


4.10591 


4.0 


2.77132 


0.0893922 


-0.582756 


0.183430 


0.0137251 


5.33827 


6.0 


3.51862 


0.0168288 


-0.269283 


0.0562110 


0.0111825 


10.1441 


8.0 


4.38982 


0.00522316 


-0.153149 


0.0238626 


0.00902933 


17.0796 


10 


2.8 














3.2 


4.57779 


0.635740 


-1.90468 


0.285430 


0.00099367 


5.15176 


4.0 


4.61313 


0.529248 


-1.51849 


0.413647 


0.00404222 


6.13667 


6.0 


4.92428 


0.112153 


-0.774269 


0.154977 


0.00651324 


10.4728 


8.0 


5.49493 


0.0339071 


-0.445769 


0.0650895 


0.00667781 


17.2497 



TABLE I: Typical values of the physical quantities of the black hole solutions in D = 4. 



Collecting all these information, we give summary of our results in the form of the relations between physical 
quantities of the black holes. Figure 3 shows various relations for neutral case (solid (black) line), charged case with 
q = 5 (dashed (red) line) and with q = 10 (dotted (blue) line). For comparison, we also give our results by the thick 
(thin) lines for dilatonic (non-dilatonic) solutions [see Eqs. (4.1)-(4.3) for the non-dilatonic case]. For the non-dilatonic 
solutions, the GB term becomes total divergence in D = 4, and the solutions are Schwarzschild and RN black hole 
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(c) 

FIG. 3: Various relations for D — 4 black hole solutions in (dilatonic) EGB systems, (a) M-th diagram, (b) M-T diagram, 
(c) M-S/Tii diagram. Solid (black) lines for q = 0, dashed (red) lines for q = 5, and dotted (blue) lines for q = 10. Dilatonic 
solutions are given by thick lines and non-dilatonic ones by thin lines. In D — 4, the non-dilatonic solutions are Schwarzschild 
(q — 0) and RN black hole (q ^ 0) solutions. 

solutions (B = 1 — 2M/r + q 2 /4r 2 ). 3 For the RN black hole, there is the minimum mass for each charge, which 
corresponds to the extreme black hole solution. There are two values of r# for the fixed charge q{^ 0) and mass (see 
Fig. 3(a)). The larger is the radius of the black hole horizon, the smaller is that of the inner horizon of the same 
solution. In contrast, in the dilatonic case, there are lower bounds on the horizon radii as well as the masses of the 
black hole solutions for given fixed charges q. They are given as (q, x, y) = (0, 2.2134, 0), (5, 2.0117, 2.7593), 
(10, 1.9686, 3.9144). These solutions are not extremal. The lower bounds on the solutions with q = and 5 are 
determined by the forbidden region in Fig. 1; the parameters of the corresponding solutions are just on the boundary 
F. The horizon radii of the smallest black holes are almost the same for this range of charge as we have pointed out 
already, while their masses depend on the charge. For q = 10, the lower bound corresponds to the singular solution 
not on F but on B2. From these figures, we see that the bigger the charge is, the further the curves move away from 
the neutral one, and the effect of the dilaton field is bigger. This is expected because the bigger the charge is, the 
bigger the dilaton coupling is. 

Figures 3 (b) and 3 (c) give relations between the mass and temperature and their entropy, respectively. For the 
neutral case, the GB term has the tendency to raise the temperature compared to the non-dilatonic solution. For 
the charged case, however, the temperature is lower than the non-dilatonic solutions. The non-dilatonic solutions 
have extremal limit with T — while all the dilatonic solutions have non-zero finite temperature. This fact gives the 
following scenario. A dilatonic black hole loses its mass by emitting radiation and continues evaporating until the 
solution reaches the minimum mass solution and the spacetime becomes singular. The entropy of the Schwarzschild 
and RN black holes in the EGB theory is different from those in GR because of the second term in the square bracket 
of Eq. (3.3). This term is proportional to a' and is absent in GR. In D = 4, the contribution from this term to 
is 1 (constant). Hence, for instance, entropy of the Schwarzschild black hole in the zero mass limit is = 1 in the 

EGB theory. 



3 The normalization of our charge is a factor 2 different from the conventional one. 
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FIG. 4: The parameter region where the black hole solution exists in D — 5. The black hole solution exists in the allowed 
region which is shadowed by (blue) thin lines, (a) The horizontal axis is x = e 1 ^ H -^' 1 ru and the vertical axis is y = e 37 **' 2 g. 
On the curve S the first derivative of the dilaton field on the horizon <j>' H diverges and the solution becomes singular. The 
curve E represents the extreme solution. The curve B gives the boundary where the solutions exist. On the curve B the second 
derivative of the dilaton field diverges at r > rn for y < 12. (b) The allowed region in terms of x and the physical temperature 
T of the black hole. The solutions on the boundary B have non-zero finite temperature. 



VI. D = 5 BLACK HOLE 



For D = 5, Eq. (2.18) to determine (j>' H reduces to 

Cr s Hl {q 2 [3C(4C - l) 7 2 + C + l] -3Cr^(3C 2 7 2 + C + l)}^ 

+ r H {(C+ l) 2 r 6 H (3Cr 6 H - q 2 ) - 2Cj 2 [9C 2 (C + l)r 12 + 3(C - 1)(2C + l)q 2 r 6 H + 2q 4 ] }^ 



+ 7 

The discriminant is 

1 



9C 2 (3C 2 + 2C- 3)r]i + 3(3C* 2 + 6C- l)q 2 r% - 2q 4 



(6.1) 



-r 2 H (q 2 - 3Cr%) 2 [l8C 6 r 12 + 30C 5 r 12 + C 4 r%(12q 2 + 5r%) + 2C 3 r 6 H (llq 2 
+ C 2 {2q 4 + 3q 2 r 6 H - 12r%) + 8Cr 6 H (q 2 + r%) + 2r]f] . 



8r 6 



H I 



(6.2) 



We show the allowed region where the black hole solutions exist in Fig. 4(a). Since the discriminant (6.2) is always 
positive for C > 0, there is no forbidden region. This is also the case for higher dimensions. There are, however, other 
bounds given by the curve S on which the first derivative of the dilaton field diverges at the horizon. There is also 
the curve E, on which the solution becomes extreme. The equations of these curves are 



S: y = 2x' 



3(x 4 + 4x 2 + 12) 
x 4 + x 2 + 48 : 



(6.3) 



E: 



V 



= 2V3x 2 , 



(6.4) 



respectively. The actual boundary of the allowed region is given by the curve B. The interval of each thin line in the 
allowed region is Ay = 0.4. The thin lines in other figures in this section have the same interval. 

For y < 12, <f>H is large near the boundary B so that the physical horizon radius Th and the charge q are very small 
as we will see soon (recall that our physical 4>h is determined by the shift symmetry (2.25)). This implies that the 
solution has the similar properties to the neutral solution. (The neural solution with higher order term of the dilaton 
field (a = 1) is summarized in the Appendix.) The neutral solution also has the lower bound on rn, where the second 
derivative of the dilaton field diverges at outer region r > rjj. We find that the same divergence occur also in the 
charged case. Krctschmann invariant diverges there. Hence the outer domain of the black hole with parameters on 
B (y < 12) is singular. For larger y (> 12.4), <j>" diverges at radius close to the horizon r w m for solutions with 
parameters on B. For larger y, the extreme curve E and S intersect at (x, y) = (2vo, 24-\/3) = (3.464, 41.57) but the 
boundary B is still below the extreme curve E. 
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(a) (b) 

FIG. 5: The parameter region where the black hole solutions exist in D = 5. (a) The diagram of the horizon radius Th and 
the charge q. For each value of the charge q, there is a lower bound for the horizon radius th- (b) The diagram of the horizon 
radius rg and the temperature T. 



In Fig. 4(b), we show the allowed region in terms of the scaled horizon radius x and the temperature T. Each 
thin curve has the same y. The allowed region is three-dimensional with another axis of y, and the boundary looks 
to turn around into inside the region. On the boundary represented by B, the temperature of the minimum size 
solution is finite, hence a small black hole may evolve to this solution through the evaporating process. For y < 12 
there may appear the singularity at non-zero radius away from the horizon. On the other hand, for y > 12.4, the 
singularity appears at the radius very close to the horizon just before the radiation stops. The qualitative difference 
of the boundary B with y < 12 or y > 12.4 mentioned above appears around x = 2.1 in the diagram. The boundary 
curve B decreases rapidly there, and no "fold" is observed for x > 2.1. 

We also give the allowed region in terms of m, q and T in Fig. 5. The curve B corresponds to the boundary 
in Fig. 4. In Fig. 5(a) we see that the solutions with the same y (depicted by thin lines) converge to the lower left 
corner corresponding to neutral solution. For example, the physical parameter of the solution with (x, y) — (2.02, 11) is 
(rn,q) = (0.177, 0.0848). The boundary B changes its qualitative property at (rjj, q) ~ (0.2, 0.1) (y 12). Figure 5(b) 
shows that the temperature gets lower in the presence of the charge but it does not vanish for all solutions. The left 
vertical boundary of B corresponds to y < 12. 

Figure 6 shows the configurations of the field functions cj>, m and S for the horizon radius r# = 2.1 for neutral and 
charged cases with q = 0, 5 and 10. This horizon radius is chosen such that the solutions exist for these charges. There 
are some different features in D = 5 dimension compared to the D = 4 case; the dilaton field decreases and the lapse 
function increases rapidly around the event horizon. Although the configurations in the neutral case are qualitatively 
different depending on the horizon radius (details can be found in Ref. [16]), we find that the variations of the field 
functions become large as the charge becomes large in general due to the dilaton coupling as in the D = 4 case. 

Here we tabulate physical quantities for the black hole solutions with the charge q = 0, 5 and 10 in Table II. 

In Fig. 7, we give the relations of physical quantities for neutral and charged solutions with different charges as 
well as those in non-dilatonic theory. For the neutral case, the non-dilatonic black hole has the zero horizon radius 
limit where M = 1 (Fig. 7(a)). It is a singular solution. There is also a lower bound on the horizon radius for the 
dilatonic solution. For the minimum solution, the second derivative of the dilaton field diverges outside the horizon. 
For the charged case, the non-dilatonic solution is described by Eqs. (4.1) and (4.3). There is the extreme solution 
for each charge, which has the minimum mass. As in the RN black holes in D = 4, the lower curve for each charge 
is the radius of the inner horizon. The graphs of the charged dilatonic solution show similar behavior of the neutral 
one except for the radius of the lowest mass solution. For the minimum solution, the second derivative of the dilaton 
field diverges just outside the horizon. 

We can see from Fig. 7(b) that the non-dilatonic solutions has zero temperature in the low mass limit regardless of 
the charge. For the charge case, the dilatonic solution behaves like non-dilatonic one. However, the temperature of 
the neutral solution raises as the mass becomes small. Figure 7(c) shows the mass-entropy relations. 
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FIG. 6: Configurations of (a) the dilaton field 0, (b) the mass function m, and (c) the lapse function S of the black hole solutions 
in D — 5. The horizon radii are m = 2.1, and the charges are q — (solid (black) line), q = 5 (dashed (red) line), q — 10 
(dotted (blue) line). 



q 


ru 


M 


Sh 


<j>H 


<t>H 


T 


S/Ei 





0.2 


0.25377 


0.512040 


7.37067 


-19.7397 


0.0563593 


0.0170506 


1.7 


2.21038 


-0.341768 


0.268320 


-1.08980 


0.0395517 


5.68795 


2.2 


3.56697 


-0.173250 


-0.168216 


-0.346471 


0.0383783 


9.84113 


4.0 


9.43189 


0.065321 


-0.305025 


0.172855 


0.0313079 


29.9771 


6.0 


19.2458 


0.0181629 


-0.164383 


0.0680353 


0.0237666 


73.5419 


5 


0.2 














1.7 


2.41969 


-0.484009 


0.678086 


-2.25502 


0.0254686 


4.86175 


2.2 


3.69719 


-0.208188 


-0.0987342 


-0.464995 


0.0352647 


9.59600 


4.0 


9.49343 


0.064932 


-0.306446 


0.171797 


0.0310818 


29.9870 


6.0 


19.2748 


0.0182388 


-0.165188 


0.0682174 


0.0237302 


73.5498 


10 


0.2 














1.7 














2.2 


4.14235 


-0.343864 


0.157902 


-1.05281 


0.0236971 


8.76096 


4.0 


9.68000 


0.0637826 


-0.310739 


0.168659 


0.0304058 


30.0171 


6.0 


19.3627 


0.0184687 


-0.167607 


0.068766 


0.0236212 


73.5735 



TABLE II: Typical values of the physical quantities of the black hole solutions in D = 5. 
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FIG. 7: Various relations for D — 5 black hole solutions in (dilatonic) EGB systems, (a) M-rn diagram, (b) M-T diagram, 
(c) M-S/Tii diagram. Solid (black) lines for q = 0, dashed (red) lines for q — 5, and dotted (blue) lines for q = 10. Dilatonic 
solutions are given by thick lines and non-dilatonic ones by thin lines. 



VII. D = 6 BLACK HOLE 



For D = 6, Eq. (2.18) to determine <j>' H reduces to 

Cr^-f^4Cr 8 H [4C 2 (6C + 7) 7 2 + C(2C + 7) + 3] - 3q 2 [4C(6C - 1) 7 2 + 2C + l] }^ 

- r H hC"i 2 [80C 2 (C 2 - C - l)r]f - 12(C - 1)(4C + l)<? 2 r^ - 9g 4 ] + (2C + l) 2 r^ [4C(C + 3)r| - 3<? 2 ] 



80C^(7C" + 2C- 2)r^ + 12(6C 2 + 6C - l)q z r% - 9q 



= 0. 



(7.1) 



The discriminant is 



-r 2 H [3q 2 ~ 4C(C + 3)r s H ] 2 [A00C 6 r]f + 480C 5 r^ 6 + 24C 4 r|(5<? 2 +6r%) + 8C 3 r s H (15q 2 - 4r%) 



C 2 (9q 4 + 36q 2 r* H + l&rfi) + 8Cr* H (3<f + Ar H H 



(7.2) 



which is always positive for C > 0. Hence again there is no forbidden region from the existence of the proper boundary 
condition in the x-y diagram in Fig. 8. There is, however, the region where solutions cease to exist if we keep y constant 
and make x smaller. The curve S where <p' H diverges and the curve E where solution becomes extreme are expressed 
as 



S: y = 2x< 



6(z 6 + 14a: 4 + 108.T 2 + 432) 



6a; 2 + 216 



(7.3) 



E: y = 2x 2 ^Q{x 2 + 2), 
respectively. The curves S and E intersect at (x, y) = (2y/5, 80V33)- 



(7.4) 
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(c) (d) 

FIG. 8: The parameter region where the black hole solution exists in D — 6. The black hole solution exists in the allowed 
region which is shadowed by (blue) thin lines, (a) The horizontal axis is x = e 7 * >ff//2 rH and the vertical axis is y = e 37 * ,ff//2 q. 
On the curve S the first derivative of the dilaton field on the horizon <j>' H diverges and the solution becomes singular. The line 
E represents the extreme solution. The curve B gives the boundary where the solutions exist. On the boundary B the second 
derivative of the dilaton field diverges just outside the horizon, (b) The magnified diagram of (a). The boundary B is separated 
from the singular curve S. (c) The allowed region in x and the physical temperature T of the black hole. There are solutions 
above the curve B on the right side. The solutions on the boundaries B have non-zero finite temperature. The boundary B 
seems to go down to the origin, (d) The diagram of the allowed region in terms of the physical quantity m and q. For each 
value of the charge q, there is a lower bound for the horizon radius rjj. 

For the neutral case, there is a solution in the zero horizon limit, while there is the lower bound for x for the charged 
solution as shown in Fig. 8(a). We find that the boundary is given by B where cf>" diverges just outside the horizon 
(r ps ru)- Kretschmann invariant also diverges there. In D = 5, the second derivative of the dilaton field diverges at 
r > rn on the boundary, but there is not such behavior in D > 6. The interval of the thin lines in the allowed region 
is Ay = 0.5. Figure 8(b) is the magnified diagram to distinguish B from S. The difference between them is small. 
Hence the boundary B is well approximated by Eq. (7.3) for x<2\/b. From Fig. 8(c), the temperature does not 
vanish on the boundary B but it is lower than that of closer solutions on the same curve. On the boundary B, the 
temperature decreases as y — > 0, which corresponds to zero charge limit, because the boundary B approaches to the 
extreme curve E in this limit. However, for the neutral solution with y = exactly, the temperature diverges in x — > 
limit. This disconnected behavior is the same as that between the Schwarzschild and the RN black hole solutions in 
GR. There are solutions above the curve B on the right side (although we do not fill in curves). The boundary B 
seems to go down to the origin. In Fig. 8(d), we display the allowed region in terms of physical quantities rjj and q. 
The left boundary in Fig. 8(d) is B, to the left of which there is no solution. 

Figure 9 shows the configurations of the field functions <fr, m and 5 for the horizon radius tr = 1-6 for neutral and 
charged cases with q = 0, 5 and 10. 

We tabulate quantitative results on the physical quantities for the black hole solutions with the charge q = 0, 5 and 
10 in Table III. 

Figure 10 show the relation between physical quantities for some fixed charges. Figure 10(a) shows the mass and 
horizon radius relation. The neutral black holes have the zero mass and zero horizon radius limit for both the dilatonic 
and non-dilatonic cases. As in the lower dimensional cases the radius of the inner horizon of the non-dilatonic solution 
is depicted. In D = 6, the horizon radius of the dilatonic solution is larger than that of the non-dilatonic solution 
for the same mass, and the difference between them is remarkable. On the other hand, the dependence of the charge 
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FIG. 9: Configurations of the functions; (a) the dilaton field <j>, (b) the mass function m and (c) the lapse function 8 of the 
black hole solutions in D = 6. The horizon radii are m = 1.6, and the charges are q — (solid (black) line), q — 5 (dashed 
(red) line), q = 10 (dotted (blue) line). 



5 


rn 


M 


Sh 


4>H 


<f>H 


T 


S/Ei 





0.2 


0.468100 


-0.0468327 


0.687743 


-0.325345 


0.209961 


0.170565 


1.7 


5.76616 


-0.205696 


0.725035 


-1.05676 


0.0492691 


14.1553 


2.2 


10.3889 


-0.182805 


0.385590 


-0.713747 


0.0454031 


29.8043 


4.0 


45.6455 


-0.00965932 


-0.157190 


0.00409228 


0.0378379 


167.850 


6.0 


129.375 


0.0230420 


-0.177248 


0.0797759 


0.0312294 


560.017 


5 


0.2 














1.7 


5.84751 


-0.215021 


0.760954 


-1.12794 


0.0472207 


13.9405 


2.2 


10.4274 


-0.184416 


0.391221 


-0.723068 


0.0450363 


29.7369 


4.0 


45.6528 


-0.00969279 


-0.157154 


0.00400513 


0.0378301 


167.848 


6.0 


129.377 


0.0230421 


-0.177253 


0.0797760 


0.0312288 


560.017 


10 


0.2 














1.7 


6.09473 


-0.253558 


0.892202 


-1.43403 


0.0402291 


13.1877 


2.2 


10.5425 


-0.189437 


0.408529 


-0.752338 


0.0439171 


29.5312 


4.0 


45.6745 


-0.00979316 


-0.157046 


0.00374374 


0.0378066 


167.842 


6.0 


129.393 


0.0230423 


-0.177268 


0.0797761 


0.0312269 


560.019 



TABLE III: Typical values of the physical quantities of the black hole solutions in D — 6. 



in the dilatonic case is rather mild. We can see that the th-M graph of each charge traces almost the same curve. 
The left end points of the curves correspond to singular solutions with the parameters on B. Figure 10(b) shows 
the relation between the mass M and the temperature T . There appear some qualitative differences depending on 
whether the solution is charged or not especially around the endpoint in the small mass limit. Figure 10(c) shows 
the relation between entropy and the mass. Compared with the rjj-M diagram, the graphs of dilatonic solutions are 
separated from each other. Here it should be noted that the entropy is not proportional to the horizon area in the 
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M 
(c) 

FIG. 10: The relations of physical quantities for black hole solutions in D = 6 EGB systems, (a) M-th diagram, (b) M-T 
diagram, (c) M-SfEi diagram. Solid (black) line for q = 0, dashed (red) line for q = 5, and dotted (blue) line for q = 10. The 
dilatonic solutions are given by thick lines, and the non-dilatonic ones by thin lines. 



higher curvature theory. The expression of entropy is given by Eq. (3.3), and the second term in the square bracket 
gives this separation. 



VIII. D = 10 BLACK HOLE 



The qualitative properties in the D = 6 to 10 cases are almost the same. However, since D = 10 is the critical 
dimension in superstring theories and important for applications, here we show the summary and some diagrams of 
the black hole solutions in D = 10. 

For D = 10, Eq. (2.18) to determine (f>' H reduces to 

Cr]fj^8Cr]f [24C 2 (70C + 19) 7 2 + 90C 2 + 57C + 7] - 7q 2 [8C(14C - 1) 7 2 + 6C + l] 

- r H {8Cr% [432C* 2 (15C 2 - C - 1) 7 2 + 540C 3 + 108C 2 + 99C* + 7] 

- 7g 2 r^ 6 [16C(C - 1)(12C + 1) 7 2 + (6C + l) 2 ] - 98CqS 2 }<P' H 

- 7 [576C 2 (99C* 2 + 6C - 4)rf 2 + 56(18C 2 + 6C - l)q 2 r]f - 49q 4 ] = 0. (8.1) 
The discriminant is 

\r% [7q 2 - 8C(15C + 7)rjf] 2 [46656C 6 r% + 34560C 5 rf + 432C 4 (7q 2 r)J i + 20rf ) 

+ 192C 3 (7q 2 r^ + 9r 32 ) + C 2 (7q 2 + 24r)|) 2 + 8C{7q 2 r^ + 12r 32 ) + 4rf ] , (8.2) 

which is always positive. Hence again there is no forbidden region in the parameter space from this restriction. 
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(c) (d) 



FIG. 11: The parameter region where the black hole solution exists in D = 10. The black hole solution exists in the allowed 
region which is shadowed by (blue) thin lines, (a) The horizontal axis is x = e l '^ H ^ru and the vertical axis is y = e Al ^ H q. The 
curve B gives the boundary where the solutions exist. On the boundary B the second derivative of the dilaton field diverges 
just outside the horizon, (b) The magnified diagram of (a). There are small differences between the boundary B, the singular 
curve S, and the extreme curve E. (c) Allowed region in x and the physical temperature T of the black hole. There are solutions 
also in the lower-right region. The solutions on the boundaries B have non-zero finite temperature. The boundary B seems to 
go down to the origin, (d) The diagram of the allowed region in terms of the physical quantity th and q. For each value of the 
charge q, there is a lower bound for the horizon radius th- 



In Fig. 11, we give the allowed region for the dilatonic black holes. The interval of the thin lines in the allowed 
region is Ay = 1. In Fig. 11(a), the extremal curve E and the singular curve S where <p' H diverges overlap with the 
boundary B, which are expressed as 



, 7(x e + lite 4 + 5712a; 2 + 164640) 
x 4 + 5bx z + 5488 



E: y = 4iV7(s 2 +30), (8.4) 

respectively. The curves S and E intersect at (x, y) = (2a/13, 162432a/574). Figure 11(b) is the magnified diagram 
of (a). We can see that there are actually small differences between these curves. On the boundary B, the second 
derivative of the dilaton field diverges just outside the horizon. Figure 11(c) shows the T-x diagram. There are 
solutions between the q = curve and the boundary B. On the boundary B the temperature is not zero. In Fig. 11(d), 
we display the allowed region in terms of x and y. Because the value of </>#(< 0.3) is small for all ranges, there is not 
much difference from Fig. 11(a). 

The configurations of the field functions (f>, m and S are qualitatively similar to D = 6 case, and we omit the figures 
here. Instead, we give quantitative results on the typical physical quantities for the black hole solutions with the 
charge q = 0, 5 and 10 in Table IV. 

In Fig. 12, we give the relations of physical quantities. There is not much difference between the dilatonic and the 
non-dilatonic solutions for these values of charge. This seems to be the tendency for higher-dimensional solutions. 
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<7 


rn 


M 


Sh 


4>H 


<t>H 


T 


S/Ei 





0.2 


0.0066645 


-0.00168778 


0.0355222 


-0.0319461 


0.996737 


0.00176109 


1.7 


6.48277 


-0.0125900 


0.163543 


-0.133113 


0.253318 


6.80565 


2.2 


276.002 


-0.0368148 


0.318372 


-0.276744 


0.125172 


593.831 


4.0 


1027.74 


-0.0497577 


0.361146 


-0.333154 


0.100354 


2787.35 


6.0 


26445.5 


-0.0693803 


0.285279 


-0.335072 


0.0646700 


115826 


5 


0.2 


— 


— 


— 


— 


— 


— 


1.7 


7.06843 


-0.0138779 


0.170767 


-0.147951 


0.221395 


6.80611 


2.2 


276.008 


-0.0368154 


0.318375 


-0.276748 


0.125170 


593.831 


4.0 


1027.74 


-0.0497577 


0.361147 


-0.333155 


0.100354 


2787.35 


6.0 


26445.5 


-0.0693803 


0.285279 


-0.335072 


0.0646700 


115826 


10 


0.2 














1.7 


8.76069 


-0.0207554 


0.193860 


-0.236732 


0.123396 


6.70390 


2.2 


276.018 


-0.0368170 


0.318383 


-0.276759 


0.125164 


593.828 


4.0 


1027.74 


-0.0497578 


0.361147 


-0.333155 


0.100354 


2787.35 


6.0 


26445.5 


-0.0693803 


0.285279 


-0.335072 


0.0646700 


115826 



TABLE IV: Typical values of the physical quantities of the black hole solutions in D = 10. 



IX. CONCLUSIONS 

In this paper we have studied asymptotically flat charged black hole solutions in the dilatonic EGB theory in 
various dimensions. The theory is the low-energy effective theory of the heterotic string. The spacetime is assumed 
to be static and spherically symmetric. The system of the field equations is so complex that it is difficult to obtain 
an analytical solution. Hence we investigate the numerical solutions. The system of the field equations has some 
symmetries which are helpful in the analysis. The results are given for D = 4,5,6 and 10, and we did not discuss 
other dimensions 7 to 9 simply because we expect that behaviors are qualitatively the same as those presented here 
from our earlier study of the system. 

We have found that there is the forbidden region on the parameter plane spanned by x and y which are the "scaled 
horizon radius" and the "scaled charge", respectively, in D = 4. Besides it, there are some boundaries (and important 
curves) of the allowed region such as the singular curve S where <p' H diverges, the extreme curve E where the solution 
becomes extreme, and B^ (i = 1, 2, 3) where the field functions and/or its derivatives diverge. The forms of the 
allowed parameter regions are different depending on the dimension, 4, 5, and 6 - 10. 

We have also studied the thermodynamical quantities. There is no extreme black hole solution with T = that 
can be obtained by taking the limit of the non-extreme solutions within the parameter range we chose. In the higher 
curvature theory, physical entropy is defined by Iyer and Wald [32]. Although there is a parameter region where 
the radius (or the area) of the black hole horizon in the dilatonic theory is smaller than that in the non-dilatonic 
theory, entropy in the dilatonic theory is always larger than that in the non-dilatonic theory. Since the dilatonic and 
non-dilatonic solutions arc the solutions in different theories, it is difficult to compare their thermodynamical stability 
and quantum transition between them in the thermodynamical sense as long as 7 ^ 0. There would be, however, 
some physical meaning in this magnitude relationship of entropy. 

It is noted again that our analysis includes the higher order term of the dilaton field which is not in our previous 
works [16-20]. To make its effects clear, we have studied the solutions for both cases a = and a = 1 in Eq. (2.4) in 
the Appendix by focusing on the neutral solutions. The qualitative properties and the relations such as M-rn, M-T 
and M-SfEi are quite similar in each case except for D = 5. In the D = 5 theory without the higher order term, the 
black hole solution with infinitesimal size exists towards rn = while there is a lower bound on the horizon radius 
th in the theory with the higher order term. For the lower bound solution, the second derivative of the dilaton field 
diverges at some radius outside the horizon, and the spacetime becomes singular. The existence of the lower bound 
of this type is confirmed also in the charged solution with y < 12 (q < 0.1). 

By our analysis it is found that the properties of the black hole solutions strongly depend on the dimension, charge, 
existence of the dilaton field. Hence both the detailed analyses of the individual systems and the investigations from 
the systematic point of view are important. 

There still remain some questions in our work. The first one is to determine precisely where the real boundary of 
the allowed region B (or B 2 ) is. It is difficult to determine it due to the fine structure of the system near the boundary 
and the numerical accuracy. The second is concerned with the fact that the temperature of the black hole solutions on 
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M 
(c) 

FIG. 12: The relations of various physical quantities for black hole solutions in D = 10 dilatonic EGB systems, (a) M-rn 
diagram, (b) M-T diagram, (c) MS /Hi diagram. Solid (black) line for q = 0, dashed (red) line for q = 5, and dotted (blue) 
line for q — 10. The dilatonic solutions are given by thick lines, and the non-dilatonic ones by thin lines. 

some boundary remains non-zero and finite. This means that the evaporation does not stop there and the black hole 
still evolves to a naked singularity or something unknown. This gives a very interesting puzzle that the singularity 
may be really formed after evaporation process and deserves further study. 

Another interesting future work is to study charged AdS black holes for application to AdS/CFT correspondence. 
We expect that this class of solutions exists for k = [17]. We leave these problems for future works. 
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Appendix A: The effects of the higher order term of the dilaton field for the neutral solutions 

In our previous work on the black holes in the dilatonic EGB theory [16-20], we examined the neutral solutions 
without the higher order derivative term of the dilaton field (a = in Eq. (2.4)). There is, however, such a term in 
general. Hence it is significant to investigate if there are any differences between the cases of a = and a = 1 . The 
neutral solutions for a = 1 can be obtained by putting q = in our model. One can find the result by looking at the 
curves with q = in Figs. 1-3 for D = 4. Comparing them with the a = case in our previous papers, we find that 
the qualitative properties and the relations such as M-rn, M-T and M-S/Yii are quite similar. This is true also in 
other dimensions except for five. Hence let us discuss D = 5 case in more detail here. 

Figure 13 displays the relations between physical quantities for neutral solutions. The black hole solutions in the 
theory with the GB term have the non-zero lower bound on their mass. Compared with the non-dilatonic BD solution, 
the mass of the lower bound in the dilatonic solutions is smaller, and the one for a = 1 is slightly smaller than that for 
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M M 



(c) (d) 

FIG. 13: Difference between a = 1 and a = for neutral black hole solutions in the D — 5 dilatonic EGB system with 7 = 1/2. 
(a) M-Th diagram, (b) M-r a diagram (The thin lines show the horizon radius), (c) M-T diagram, (d) MS /Hi diagram. The 
solid (red) line is for o = 1, the dashed (blue) line for a — 0, the dotted (green) line for Tangherini solution (higher-dimensional 
generalization of the Schwarzschild solution without the GB term), and the dotted-dashed (black) line for the BD solution (in 
non-dilatonic EGB theory). 

a = (Fig. 13(a)). We also find that for the a = case the black hole solution with infinitesimal size exists towards 
Tn = while there is a lower bound on the horizon radius r# for a = 1. This bound corresponds to the boundary 
B in Figs. 4(a) and 5(a). With the parameters around the lower bound B the second derivative of the dilaton field 
grows suddenly at some radius r > ru- This behavior is shown in Fig. 14 for the parameters just before touching the 
boundary. It is expected that <j)" should diverge on the boundary. Since the Kretschmann invariant also diverges at 
that radius, the solution of the lower bound is singular. If we compare black holes with the same mass, the horizon 
radius for a = 1 is a little larger than that for a = 0. As a general tendency, values of the dilaton field <fi itself and its 
variation <\>' are larger in magnitude for a = 1 than those for a = 0. For the large black holes (rjj > 1), the differences 
are indistinguishable. 

1000 | 1 , , 1 1 

800 - 
600 - 
400 - 




0.5 1 1.5 2 

r 

FIG. 14: Configuration of <f>" for parameters q = 0, rn ~ 0.103555 just before touching the lower bound in D — 5. The peak 
is at r = 1.702. This behavior is basically the same for the charged solutions with y < 12 (q < 0.1). 



Figure 13(b) shows the radius of the curvature singularity (thick curves) inside the black hole horizon (of which 
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radius is depicted by thin curves). For M < 7.5, the curvature singularity is located at the center of the black hole. 
For the large mass M > 7.5, however, the singularity is expanded to a finite radius r s > 0. We called it a fat 
singularity [20]. Even if the higher order term of the dilaton field is included, this tendency is unchanged. 

The temperature T is slightly higher in a = 1 than in a = case for M > 0.6, while the order changes near the 
minimum mass (Fig. 13(c)). Although it might appear that the difference is very small in D = 5, the difference is 
further smaller in other dimensions. In addition, the temperature of the minimum size solution is finite, hence a 
small black hole for a = 1 may evolve to this solution through the evaporating process, and the naked singularity 
might appear at non-zero radius. Such behavior cannot be seen for higher dimensions even if the higher order term 
is included. Figure 13(d) is the entropy vs. mass diagram. 
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